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In the study of d-dimensional quantum channels (d > 2), an assumption which is not very restric- 
tive, and which has a natural physical interpretation, is that the corresponding Kraus operators 
form a representation of a Lie algebra. Physically, this is a symmetry algebra for the interaction 
Hamiltonian. This paper begins a systematic study of channels denned by representations; the 
famous Werner-Holevo channel is one element of this infinite class. We show that the channel 
derived from the defining representation of SU (n) is a depolarizing channel for all n, but for most 
other representations this is not the case. Since the Bloch sphere is not appropriate here, we 
develop technology which is a generalization of Bloch's technique. Our method works by repre- 
senting the density matrix as a polynomial in symmetrized products of Lie algebra generators, 
with coefficients that are symmetric tensors. Using these tensor methods we prove eleven theo- 
rems, derive many explicit formulas and show other interesting properties of quantum channels in 
various dimensions, with various Lie symmetry algebras. We also derive numerical estimates on 
the size of a generalized "Bloch sphere" for certain channels. There remain many open questions 
which are indicated at various points through the paper. 
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I. INTRODUCTION 

It has been known for more than three decades through the work of Kraus (1971), and that of Choi (1975), that 
a general channel admits an operator-sum representation, so its action on an n x n density matrix p has the form 
J2 M M^pM^. A special case occurs when the M M are Hermitian and maximal in number, so that p — . . . n 2 — 1. 

Define Ao = 1 and let A a for a > 1 denote the n-dimcnsional Gell-mann matrices, which are standard generators 
for the Lie algebra su„. Then {A M } is a basis for the space of Hermitian matrices. Taking M proportional to the 
identity, there exist constants U a b such that 

M a = Uabh- 
b 

If U is a unitary matrix, then we may take each M M proportional to A^ without changing the quantum channel defined 
by these Kraus operators. 

In the latter case, one may readily calculate M^pM^, since p itself may be expanded in the A M basis, and the 
Gell-mann matrices in any dimension satisfy elegant product identities. Indeed, we calculate this operator explicitly 
in Sec. IV. Exploiting the product identities leads to a simple, elegant and explicit formula for the action of the 
5U„ channel; if p — ^1 + X^ 1 V'V tnen tne cnannc l multiplies v by a scalar. Thus the assumption that the Kraus 
operators are proportional to generators of su„ is extremely useful as a calculational tool, but it is not merely that. 
A quantum channel models the interaction of a decohering system with its environment, and the identification of the 
M M as generators of a Lie algebra is related to a symmetry of the interaction Hamiltonian. See (Lidar et al., 1998) 
and (Lidar and Whaley, 2003) for details. 

The generalizations of the calculations in su„ described above to other Lie algebras and to higher-dimensional 
representations are illuminating, and have not appeared in the literature before. These cases necessarily have the 
property that not all of the Hermitian matrices in that dimension are linear combinations of representation matrices, 
so direct generalization of the calculational method outlined above for su„ will not work, and a new idea is required. 
This is the subject of Section V. We give a particularly detailed analysis of the three-dimensional representation of 
SU2; however, many of the formulas we use there generalize readily to higher spin. In an interesting twist, the spin-one 
case turns out to be a generalization of the Werner-Holevo channel, and suggests that Lie algebra channels may play 
an interesting role as examples or counterexamples to the well-known AHW conjecture. 

Section VI analyzes the channel based on the exceptional algebra 02- In Section VIII we generalize some aspects 
of the Bloch sphere to density matrices constructed from Lie algebra representations. It is shown that for each 
representation, there is a class of density matrices parameterized by a closed, bounded (hence compact) submanifold 
of Euclidean space, which we term 'the Bloch manifold.' Explicit bounds are given on the size of these manifolds. A 
general method is given for finding the Bloch manifold exactly, using trace identities. 

II. KRAUS DECOMPOSITION 

A quantum channel is a model for a particular snapshot of the time evolution of a density matrix, and especially 
for the evolution of pure into mixed states. Let TL be a Hilbert space of dimension n, and let gl(TL) denote the vector 
space of all linear transformations from TL — > TL. A map £ : gl(TL) — ► Ql(TL) is said to be completely positive ( CP) if it 
is linear and £ ® 1 is positive on TL ® TL 1 for all TL' . The map is said to be trace-preserving if Tr £ (A) = Tr A for all 
A&q[{TL). 

Definition 1. A CPT map or stochastic map is a completely positive, trace-preserving linear transformation £ : 
S l(TL)^ 5 l(TL). 

Of central importance to the current work is a famous theorem of Kraus (1971) which gives a simple representation 
of any quantum channel. 

Theorem 1 (Kraus Decomposition Theorem). For any quantum channel £, there exists a finite set of operators 

M ,M 1 ,M 2 ,...,M k , where k < (dimTL) 2 , 

such that 

£(p) = Y i M liP Mt with Y. M t M » = 1 - (!) 

n ft 

In this situation, (1) is called the Kraus representation, the operator sum representation or the Stinespring form ; 
while ^^M^M^ = 1 is sometimes called the normalization condition and is just the statement that the map is 
trace-pres erving. 
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A proof of this theorem may be found in the original article of Kraus (1971), or in the book by Nielsen and 
Chuang (2000). We simply note that the converse, namely that any operator of the form (1) satisfies the conditions of 
Definition 1, is clearly true. A stochastic map may also be obtained as the partial trace of a unitary conjugation on a 
larger space; see (Ruskai, 2002, Sec. III.D) for a discussion. The representation (1) is sometimes called the Stinespring 
form since its existence follows from the Stinespring dilation theorem (Stinespring, 1955). 

This is a general framework, and in order to obtain explicit results, further assumptions are necessary. A mathe- 
matically elegant assumption is that the possible errors introduced in the decoherence process are not arbitrary, but 
that they correspond to the action of the infinitesimal generators of a Lie group G of continuous symmetries. This 
provides a simple model for symmetry breaking in quantum mechanics. 

The situation just described, in which the error generators are also generators for a matrix representation of a 
semisimple Lie algebra, follows naturally from the model of Markovian dynamics considered by Lidar, Chuang, and 
Whaley (1998). This was shown to have important consequences for the possibility of decoherence-free dynamics; see 
(Lidar and Whaley, 2003, and references therein) for an up-to-date review. The present work may be considered as 
a further exploration of the consequences of that model, for a snapshot of the time evolution. 

The qubit depolarizing channel is a model of a decohering qubit in which the decoherence has an SU(2) symmetry. 
With probability 1 — p the qubit remains intact, while with probability p an error occurs. The error can be one 
of three types, each equally likely. These errors are implemented by applying Pauli matrices to the qubit state. In 
other words, an error involves applying one of the generators of the Lie algebra SU2 to a vector in its irreducible 
two-dimensional representation. These generators admit direct physical interpretations as bit-flip errors, phase-flip 
errors, or combinations of those. 

The qubit depolarizing channel admits a generalization to a channel with k possible errors based on an n-dimensional 
representation H of a fe-dimensional Lie algebra g, for which detailed properties have not been previously investigated, 
and which is the main topic of the present work. As we develop the general theory of these channels in the following 
sections, we will see that not all Lie algebras can give quantum channels (at least not in the way outlined here), and 
even for Lie algebras which do give channels, not all representations are acceptable. For g semisimple, it is necessary 
that the quadratic Casimir take a single value on all elements of the representation space. This holds for all irreducible 
representations, and some reducible ones. The non-semisimple case is more difficult, and its treatment will be deferred 
to a separate paper. 

III. QUANTUM CHANNELS FROM LIE ALGEBRA REPRESENTATIONS 

This section contains our notations and conventions for the generalized depolarizing channels which will be studied 
in detail in later sections. The possibility of defining a quantum channel based on a representation of a compact Lie 
algebra was mentioned briefly, but never elaborated upon, in a paper of Gregoratti and Werner (2003). In any case, 
it is not necessary that the Lie algebra be compact. 

A. Pure Lie Algebra Channels 

It is a standard convention (Georgi, 1982; MacFarlane et al., 1968) to normalize the canonical generators for the 
defining representation of su„ so that 

Tr (A a A 6 ) = 2S ab . (2) 

This has the desirable feature that the canonical generators for n — 2 are the Pauli matrices, and those for n = 3 arc 
the familiar Gell-mann matrices, while inserting factors of 2 in certain formulae. With convention (2), these generators 
will be orthogonal but not orthonormal with respect to the Killing form. We return to this point below. 
On a general semisimple Lie algebra, the Killing form K is defined as 

K(X, Y) = Tr(ad(A) o ad(F) : q - g), 

where the trace is taken in the adjoint representation. At the moment we focus on semisimple algebras g, for which the 
Killing form is nondegenerate, and return to treat non-semisimple algebras in a later section. Let a be an irreducible 
representation of g, let Xi be any basis of g, and let X[ denote the dual basis with respect to the Killing form. The 
Casimir operator 



tf 2 (a)=5>(*i)a(*;) 

i 



4 



does not depend on the choice of basis, and by Schur's lemma is proportional to the identity, so we write 
C-2,{a) = c 2 (a)l. If Xi is orthonormal with respect to K, then C 2 (a) = J2 i a(X i ) 2 . For reducible representations, 
C 2 (a) may not be proportional to the identity. 

Definition 2. Let q denote a Lie algebra of dimension k, with basis {Xi : i = 1, . . . , k}. Let a be an irreducible 
^-representation on the Hilbert space 7i. The generalized depolarizing channel or Lie algebra channel is defined to be 
the channel in which an error occurs conditionally with probability p, causing an initial state £ 7i to evolve into 
an ensemble of the k states a{Xi) all with equal likelihood. 

The Kraus operators for the channel of Definition 2 are given by 

M = ^Y~p 1, Mi = ^Xp a(Xi) . (3) 

where A is a normalization constant which will be fixed momentarily. The operators M M are hermitian if the rep- 
resentation is unitary and if p e [0, 1], and are constrained to satisfy M M M M = 1, which fixes the value of the 
constant A appearing in (3). By definition, 

X;^ = (i-p)i+Apj;a(jr j ) 3 . 

fj. i 

If J2i a (Xi) 2 — Z ■ 1, where Z is a constant (which in most cases we can take to be real), then 

If Xi is orthonormal with respect to the Killing form, then Z = c 2 (a) . More generally, if the basis satisfies 

K(Xi,Xj) = nSij , n > 0, 
then it can be rescaled to an orthonormal basis by a single constant. In this situation, 

Z = nc 2 (a), A=i. (4) 

Defining the Killing norm by ||x||^ = K(x, x), we note that if 

IPQIk^U^lk, 

for some pair of indices i,j, then the normalization condition cannot be satisfied. 

What if the representation is reducible? Suppose H = V W as a direct sum of irreducible g-modules, and Xi is 
orthonormal with respect to K. Then there exist independent constants Zy and Zw such that the operator 

c a («) = E«(*o a = (t I) 

i v ' 

as a block decomposition on V W. If Zy 7^ Zw, then it is not possible for the Kraus operators (3) to give a trace- 
preserving map. On the other hand, if Z v = Z w then they do define a CPT map even though the representation is 
reducible. 

What if p > 1? Then M = i^/p — 1 1, and we have 

£M„Mt = (2p-l)l. 

fj- 

Thus the map cannot be trace-preserving unless p = 1, which is a contradiction. A similar argument shows that p < 
does not give a trace-preserving map. Thus, if we wish to study the framework of Definition 2, then we must limit 
ourselves tope [0,1]. 

We summarize the results of the last few paragraphs in a Theorem. 

Theorem 2 (Normalization). Consider the Kraus operators 

M = (1 - p) 1/2 l, and Mi = (Ap) 1 / 2 ^), 



for i = 1 . . . k . If 
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(i)p& [0,1], 

(nj TTie representation a of g is a direct sum of irreducible representations all with the same quadratic Casimir, and 
(Hi) 3 n > such that K(Xi,Xj) = n5ij for all i,j, 

then M^M^ — 1 with A given by eqn. (4). Conversely, if any of (i)-(iii) is not satisfied, then (except in trivial 
cases) there does not exist A s.t. J^p. M^M^ = 1, and the M's do not give rise to a quantum channel. 

The coefficients of the in (3) admit a natural "probability of error" interpretation, but in Section III.B wc 
investigate the possibility of modifying them to complex coefficients in order to obtain a new channel. We find that 
no new channels arise unless one is willing to promote the coefficients to operators. 

Using (1), the Lie algebra channel has the explicit Kraus decomposition 

k 

p^£(p) = (l-p)p+|^«(A J )pa(A,). (5) 

i=i 

As is proven in standard textbooks (Georgi, 1982, see Theorem 8.9), the trace of any generator of any representation 
of a compact simple Lie algebra is zero, so in particular, the a(Xi) are traceless. Moreover, it is clear that this 
transformation satisfies the defining properties for a quantum channel, given here as Definition 1. 
Two operator-sum representations 

Y^M^pMl and £ N v pNt 

describe the same channel if and only if there exists a unitary matrix U uli such that N v = U^^M^. Therefore, it 
is immaterial which basis of the Lie algebra that we use, as long as the two bases are related by a U(N) similarity 
transformation. As noted in Theorem 2, in order to build a channel satisfying the normalization condition, we are 
forced to use a basis satisfying "orthonormality," K(Xi,Xj) = ndij. But any two "orthonormal" bases in this sense 
are related by a unitary transformation, so the CPT map constructed above is independent of the basis chosen for g. 

Given a Lie algebra g and a representation a on a vector space of dimension d, the CPT map (5) is a model for 
decoherence through a cZ-level noisy quantum channel, with errors that are not completely arbitrary; rather, they 
transform the state in a way determined by the representation of g. 

The channels (5) have an extremely interesting structure. For a certain subclass of possible Lie algebra representa- 
tions, the channel (5) has an action which, like the qubit case, is most simply described by a Bloch parameterization 
with polarization vector v G R fc , where k = dimg. In these cases, we show that (5) decreases the length of v, and so 
deserves the title 'generalized depolarizing channel.' In other cases of interest, a single Bloch vector is not sufficient, 
but the action of the channel can be described by similar rescalings of symmetric 2-tensors or higher-rank objects. 

A natural step, which we begin in the next section, is to calculate the expression (5) explicitly in certain represen- 
tations of classical Lie algebras. 

Remark 1. When we use the terminology "the g-channel," where g is a semisimple Lie algebra, the fundamental repre- 
sentation of g is implied. Examples of fundamental representations include the n-dimensional defining representation 
of su„, and the 7-dimensional irrep of G*2. 

It is easy to see that the Lie algebra channel (5) always has the property of being doubly stochastic, i.e. £ (1) = 1. 
See for example (Gregoratti and Werner, 2003) for further discussion. 

B. A Note on Coefficients and Extensions 

As discussed prior to Theorem 2, for p £ [0, 1] the channel defined by (5) is CP but not T, and it is possible to 
recover a CP channel only if we consider different coefficients for the Kraus operators (3). To this end, let us first 
consider 

Tfl 

M Q = m l and Mj = — a(Xi) , (6) 

v Z 

where mo, fh G C are some complex constants. Then to obtain a trace-preserving map, we require 

^TMtM M Hmo| 2 + H 2 = l. 
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This condition is equivalent to the statement that the point (mo, to) e C 2 = R 4 lies in the unit 3-sphere S 3 cl 4 . 

We can now view the coefficients of the Kraus operators (3) as the projection S 3 — ► S 1 . Introduce a parameter 
q E [—1, 1] such that p = q 2 , and write (3) as M = ± \/l — q 2 1, and Mj = (q/y r Z)a(X i ). Then ignoring y/Z, the 
coefficients of Mi and Mo give a point on the unit circle. Further, too and fh only enter through the square of their 
magnitude, so the two additional parameters associated to projecting from the 3-sphere are fictitious, and (5) is in 
fact the most general channel of this kind. 

A non-trivial generalization is obtained by promoting too and fh to operators. However this "generalization" is a 
special case of a well-known operation which extends an existing channel £g using any set of operators which satisfy 
the normalization condition (1). Given two sets of Kraus operators A\, . . . , A r and B\, . . . , B s acting on the same 
vector space and satisfying 

r s 

»=i j=i 

we note that the set of operators 

{A 1 ,...,A r - 1 ,B 1 A r ,...,B s A r } (7) 
also satisfies the normalization condition, because 

s 

^ ] ( Bj A r ) ^ Bj A r — Aj r A r . 
i=i 

This construction is natural with respect to the channel £b defined by in the sense that if {S-} is another set of 
Kraus operators defining the same channel, then the channel defined by (7) is also the same. Naturality does not hold 
for the A operators, but this will not concern us here. We call this procedure the extension of £b by the A-operators, 
on the element A r . 

For example, one may notice that the operators Z~ x l 2 o.{Xi) of the previous section satisfy the normalization 
condition since the sum of their squares is a Casimir element, and the normalization constant Z cancels the numerical 
factor. Consider this the B-channel, and extend it on every element by the same set of Kraus operators. This yields 
a "double g-channel" with Kraus operators 

{^XiXj : i,j = l...k}. (8) 

These operators generate the image of under the universal homomorphism expressed in the commutative diagram 
( 18 >\ 

This underscores the fact that, aside from the basic examples of new quantum channels provided by Section III. A, 
many further examples may be obtained by extension, as in (7). As in the basic Lie algebra channels, computations 
with extended channels are facilitated by the existence of non-trivial identities which exist among the representation 
matrices. Channel (8) is interesting because for many representations, the matrices a(Xi) do not span the entire 
space of traceless dx d matrices, but the set of products a(Xi)a(Xj) spans a subspace of larger dimension. Therefore 
the extension leading to (8) is a way of generating a channel whose Kraus operators come closer to spanning the 
space of all matrices in the appropriate dimension. If a density matrix were written as p = Y]jj WijXiXj, and if the 
representation satisfies an identity for reduction of products of six generators, then one can calculate the action of (8) 
on p explicitly. 

We are now in a position to interpret the channel defined by (6) with complex coefficients mo, fh as the extension (7) 
of the nontrivial Lie algebra channel Bi = Z^ x / 2 a{Xi) by the identity channel with the unusual Kraus representation 
A\ = m l, A 2 = ml. In case fh = q e [—1,1] and to = ±\/l — q 2 we recover (5). Given any channel whose set of 
Kraus operators do not contain 1, we can always extend it so that they do contain the identity by this method. 

For the rest of this paper, we will assume that the Kraus operators take the form (3) in order to retain the beautiful 
probabilistic interpretation given by Definition 2. As we continue, we will keep the fact in mind that extensions are 
possible, and develop methods which easily generalize. 

IV. THE SU„ CHANNEL 

The su„ channel, our first example, is the channel built from the n-dimensional defining representation (also called 
'standard representation') of su„. It is simpler than most other channels studied in this paper, because it admits a 
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complete solution. Its action on any arbitrary input density matrix can be calculated in closed form using the Bloch 
parameterization, and in all cases it is a depolarizing channel. 

One reason for the beauty and simplicity of the su n channel is that any n-dimensional density matrix admits a 
Bloch vector parameterization in terms of su„ generators. This is because k = dim(su„) = n 2 — 1 is only one less 
than n 2 , the dimension over M of the space ofnxn Hcrmitian matrices. 

Any n x n Hermitian matrix p may be represented as 

p= i(tr(p)l + T), Tesu n , 
and having chosen a basis X a for su„, it follows that 

k 

T = ^v a X a = vX, 

a=l 

for some coefficient vector v. In analogy with the well-known parameterization of the 2x2 density matrices as the 
interior of a sphere, we will refer to v as the Bloch vector. 

For n > 3 it may be hard to visualize the geometry of the space of density matrices in terms of the geometry 
of v. This question was first considered in the n — 3 case by MacFarlane et al. (1968). Section VIII undertakes a 
systematic study of the geometry of the space of v which lead to a valid density matrix in various representations. 
We call this space the Bloch manifold and give details of the geometry for a number of important examples, including 
all representations of SU2, and the n-dimensional irrep of su„. 

In this section, we take a to be the standard representation of su„ on a vector space 7i of dimension n. For simplicity, 
we let Xi denote both the generator of su„ and its image under this representation. One could now compute the 
quadratic Casimir in the standard way using roots and weights, but it will turn out that the value of this Casimir 
as well as all other properties we will need to obtain a complete solution to the su„ channel follow from the single 
relation 

XiXj = f36ij 1 + ^ QijkXk ■ (9) 
k 

for some constant [3 and tensor Qijk- Of course, this relation is just the decomposition of a Hermitian matrix into a 
trace part with trace nf36ij, and a linear combination of the Xk, which generate the space of traceless matrices. 
Elements of the standard basis of su„ are called Gell-mann matrices, and they satisfy 

TriXiXj) = 25ij, 

so (3 = 2/ n. Many properties of the Q tensor already follow from the single assumption that Xi generate a Lie algebra. 
It is immediate that Q[ij]k = ifijk where [ij] denotes antisymmetrization, and fyk is 1/2 times the structural tensor 
of the Lie algebra. It follows that 

Qijk — d{jk ~~\~ ifijk 

for some dijk symmetric in the first two indices. Also, (9) implies 

{Xi,Xj} = -S i: jl + dijiXi. 
n i 

Multiplying by X k and taking the trace yields 

dijk = j Tr({Aj, Xj}X k ), 

therefore the d-tensor is completely symmetric, and interchange of any two indices has the effect of complex conjugating 
Q. Since ^ XiX t is a multiple of the identity, 

dak = \ Tr X t X t )X k ) = (10) 

i i 

It follows from the associativity of matrix multiplication that 

2 
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with a sum over m implied. Contracting j and k and using (10) yields 

4 

dijmdljm — fijmfjlm ~l~ (2il , 

By a general property of compact semi-simple Lie algebras, the structure constants satisfy 

fijkfijk = n5 l i , (11) 

4 

Therefore, dij m dij m = (n )Su. Using this and (11), we obtain 

4 

For this basis of su„, Z = 2k /n, where k — n 2 — 1. The action of the channel 

k 



2k 



on the density matrix 



is given by 



p=Utr(p)\ + v-X) 



£{p) = + • X + !L VjXtXjXt . (13) 



2k 

1,3 



Using (9) to expand the triple product, we have 



Y VjXiXjXi = (3v ■ X + VaQijil + Y VjQijkQkiaXa . 

j\i i,3 i,j,k,a 

Since £(p) has unit trace, it must be the case that Qiji = 0. The same conclusion also follows from (10), but it is 
amusing to see that ^ ■ Qiji must vanish because this is a CPT map. Therefore, 



£{p) = + 1 P + P/k v .1 + 4^ Vj Qi jk Q kia X a . (14) 

Using (12), we have finally 



n 2k 

i,j,k,a 



£( Pv ) = -(tv(p)l + f(p,n)vX), 
n 



where 



p (1 - p)n 2 - 1 

f(p,n) = l-p--= n2 _ i . (15) 



In the qubit case, f(p, 2) = 1 — 4p/3, which is consistent with standard results. 

The su„ channel maps an initial density matrix to a linear combination of itself and the identity, i.e. it has the 
form 

A A (p) = Ap+(i-^)l. (16) 

This is the standard definition of the n-dimensional depolarizing channel. The information-carrying capacity of 
this channel was studied in great detail by King (2003), where notably the Amosov-Holevo- Werner conjecture was 
established for channels which are products of a depolarizing channel with an arbitrary channel. Channels based on 
representations of semisimple algebras generically do not take the form (16), except possibly on special subsets of 
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the space of density matrices. See Section V and in particular Theorem 3 for a Lie algebra channel that is not a 
depolarizing channel. 

The depolarizing channel on an ra-dimcnsional Hilbert space satisfies complete positivity if and only if 

! - < A < 1 . 



1 

The su„ channel has the form (16) for A = f(p, n). Note that the relation 

-1 



n 2 - 1 



</(p,n)<l 



holds for all n>2. In fact, f(p, n) saturates both of these inequalities at the endpoints of the allowed range, < p < 1. 
At the special value p=l — n~ 2 , the siin channel is a constant map from K" _1 into the space of density matrices: 

£{p v ) = —1 for all v, at p = p c = 1 — n~ 2 . (17) 

Physically, if the probability of error happens to be p = p c , then su„-decoherence evolves an arbitrary initial density 
matrix into a completely uniform ensemble consisting of pure states with equal probabilities. This is the "worst" value 
of p, in the sense that all information about the initial density matrix has been lost. This result is stable in the sense 
that if p is only approximately equal to the critical value, the initial density matrix decoheres into an approximately 
uniform ensemble. 

We will see in Sec. V that for other Lie algebra channels, there are multiple critical values of p which generalize 
(17); this analysis leads to an interesting decomposition of the space of density matrices on H which is discussed in 
Section V.F. 



V. OTHER REPRESENTATIONS 
A. General Remarks 

In the n-dimensional standard representation of su„, the representation matrices a(Xi) span the space of all 
traceless Hermitian matrices, and thus an arbitrary initial density matrix can be expressed in terms of the a(Xi) 
and the identity. As we consider higher-dimensional representations, the representation matrices become increasingly 
sparse in the space of all traceless matrices, and thus only some fraction of the set of all possible density matrices 
can be expressed in the form d~ l l + ^2iVia(Xi). This is not all the bad news; for higher-dimensional irreducible 
representations (irreps), there is generally no analogue of the identity (9) which holds for su„. 

Therefore, the simple calculations we have done for the n-dimensional irrep of su„ do not generalize in any simple 
way to other representations; new ideas are needed. In this section, we develop methods for dealing with the general 
case of arbitrary representations. Let d — d a denote the dimension of the representation a, and gl d as usual denotes 
the associative algebra of all d x d matrices. 

A representation of g lifts to a unique associative algebra homomorphism <f> of the universal enveloping algebra 
U(q), by the universal property most elegantly expressed in the commutative diagram 




The action of <fi is simply to convert the tensor product to matrix multiplication, i.e. <f>(x ® y) = 4>{x) ■ <p(y), etc. The 
interesting property about this commutative diagram, and one which gives a computational method for Lie algebra 
channels, is that if <j> is an irreducible faithful representation and if jj is a semisimple Lie algebra, then <f> is surjective. 

This surjectivity has the consequence that for any representation of said Lie algebra, every density matrix can be 
represented as a linear combination of products of the representation matrices. In other words, the new calculational 
method outlined in this section will always work. Before continuing our discussion of this, let us consider a simple 
but nontrivial example, the spin-1 channel, in complete detail. 
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B. The Spin-1 Channel 

Consider the spin-1 representation of su 2 . We use standard angular momentum notation, in which 

, /0 10\ 1 /0 -i 0\ (l 

Ji = -= 1 1 , J 2 = -= \i -i , J 3 = 
v 2 \0 1 0/ v 2 \0 i J \0 -1, 

Before generalizing to arbitrary density matrices, we restrict attention to the simpler example of density matrices p 
which arc of the form 

p v = ^{l + v-J), veR 3 . (19) 



Then 

fJn.A = ll 4- - 

3 6 



£( Pv ) = ll + ^—E v .J + ^J2 v bJ a JbJa- (20) 



a, b 



The relation analogous to (9) does not hold, i.e. J a Jb is not a linear combination of 1 and { Ji : i = 1 . . . 3}. In this 
special case, the triple product appearing in (20) simplifies considerably; 

JaJbJa = SabJa (nO Sum), (21) 

which implies that 

£( Pv ) = l(l+(l-E)v. J) . (22) 

This takes the form (19) with v — > (l — | )v. Thus, for 3 x 3 density matrices admitting a Bloch parameterization, if 
p is a probability then the spin-1 channel scales the Bloch vector by a number between 1/2 and 1. 

Interestingly, we can go further and find a Bloch-type picture of the spin-1 channel on a general density matrix. 
The six elements of the form 

J(aJb) = ^ {JaJb + JbJa) , 

together with Ji, J 2 , and J 3 , span the space of 3 x 3 matrices. Therefore an arbitrary 3x3 density matrix p can be 
written as 

Pv,w = v ■ J + ^2wabJ( a Jb) (23) 

a. b 

for some vector v and symmetric tensor w. 

We use standard physics normalizations which entail that for the spin s representation in d = 2s + 1 dimensional 
space, 

Y, J a= X1 ^{JaJb) = ^ 5 ab . 

a 

where A = s(s + 1). Then we have 

It follows that in order to have a density matrix, we require tv(w) — 3(dA) _1 . For s = 1, tr(u>) = 1/2. 

Theorem 3 (Action of the Spin-1 Channel). The action of the spin-one channel on the vector and symmetric 
tensor are v — > v' and w — > w' , where 

Va - V a ' = (1 - |K, (24) 
W a6 -» W a fc = (1 - yj^ah + -Oab • 
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Proof. The asymmetric quadruple product identity 

^^JiJjJkJi = Sj k J 2 — J k Jj (25) 

i 

implies the symmetrized identity 

^ JiJ{ a Jb)Ji = Sabf 2 - J( a Jb) ■ (26) 

i 

Using the latter and (21), a straightforward calculation shows that £(p v , w ) is equal to 

(l - |)u • J + (l - y) ^w ab J {a J b) + ptv{w)\ , 

a.b 

which implies the stated result, since for spin-1, we have tr(w) = 1/2 and 1 = J2 a b \^abJ( a Jb)- d 

We refer to identities of the form (25) as "4 — > 2 identities," because they relate degree 4 polynomials in the 
generators to degree 2 polynomials. We have also seen one "3 — > 1 identity" in equation (21). 

It is possible to iterate formula (24), with interesting results. Clearly, after n applications of the channel, v — ► 
(1 — | ) n v. Consider a W-state, i.e. a state of the form 

Pw — ^ WabJ(aJb) , 
a, b 

and let £ n denote n applications of the spin-1 channel. 

Theorem 4 (Iteration Formula) . The action of £ n on w is the following: 

w-pr^ F {n) (p){l - 6w) + w 

where F^(p) is a degree n polynomial in p, determined as follows. F^(p) = 1 — 3p/2, and the for n > 1 are 
determined by the recursion relation 

F(-+ 1 )(p) = (l-|)jrW (p ) + £. 
Interestingly, this recursion relation has the same coefficients as the transformation (24) of w itself. 



C. Pure States of the Spin-1 Channel 

In any number of dimensions, one can find a class of pure states in the following way. Let a £ R™, and consider the 
symmetric n x n matrix = OjOj. Then evidently, 

P 2 =a 2 P, and Tr(P) = a 2 . 

It follows that 

P 2 = P & Tr(P) = 1 aeS n , 

where S n denotes the n-dimensional sphere. 

If P is a density matrix, then it is a pure state. With n = 3, these pure states are precisely the pure states that 
arise from the symmetric term in (23), assuming we take the most convenient choice of basis; i.e. the one in which 
the generators for the spin-1 representation are 

/0 \ / i\ /0 -i 0\ 

Si = -i , S 2 = , Ss = \i . 
\0 i / \-i 0/ \0 0/ 



12 



Since density matrices of the form 

Pw = ^ w a bS( a S b ), 

a,b 

arise in the Bloch-type parameterization for spin-1, it is natural to ask when this type of density matrix is pure. 
Solving the equation 

(p w ) 2 = Pw (27) 

for the components of w, we find several two-parameter families of solutions, and a one-parameter family of solutions. 
For the two-parameter families, up to signs we have p w = P, with the a-vector given by 



ai = \/w 2 2 + w 33 , a 2 = \J \ - w 2 2 , a 3 = y I - W33 . 

The off-diagonal coefficients of w are also determined in terms of W22 and w 33 , so there are indeed only two free 
parameters. 

More precisely, p w = P is one family of solutions; the others are obtained by changing the signs of any two of the 
off-diagonal components of P above the diagonal, and changing corresponding signs below the diagonal so that P 
remains symmetric. Hermiticity of P requires that a must be a real vector. This means that 

1 1 

w 2 2 < 2 ) w 33 < - , w 2 2 + W33 > , 

which constrains the point (w22,w 33 ) to lie in the interior of a certain triangle; see Figure 1. 



0.75 




-0.75 



Figure 1. 

In addition to the two-parameter families discussed above, there is a one-parameter family of solutions to (27) 
taking the form 



/ I 

2 









\2 







1 



— UJ 



where u — w 33 . 



D. Higher Spin Representations 

Note that the triple product (21) and quadruple product (25) identities are simply certain elements of the ideal 
= ker(</5), where <f> is the representation of the universal enveloping algebra, as in (18). The larger this ideal, the 

more product identities there will be in the representation of interest. For higher spin, we have the following 3 — > 1 

identity in the spin s representation of su 2 , 

3 

53 J i J aJt = (A - 1) J a , where A = s(s + 1) . (28) 

i=l 
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There is also a generalization of the 4^2 identity valid for general spin s, 

3 

JiJaJbJi = (A — 2)J a Jb + XS a bl — JbJa ■ 

i=l 

The latter has the more convenient symmetrized form: 

3 

^ JiJ( a Jb)Ji = (A - 'S)J(aJb) + A<5 afc l . (29) 

i=l 

Theorem 5 (Higher Spin Channel). Let Ji, J2, J3 be canonical generators for the spins representation of SU2 in 
dimension d = 2s + 1, and let £ s denote the spins channel. Defining p ViW = v ■ J + J2 a b w a bJ( a Jb), we have 

3 

to{p v ,w) = 1 tr(io) = — 
where A = s(s + 1). The action of the spins channel is 

£ s (pv,w) = - ^jv ■ J + (l - ^pjw a bJ( a Jb) +ptr(io)l . 

Proof. A straightforward application of (28) and (29). □ 

It is now clear that the action of the spin channel is more complicated than the scaling of a single vector. It is 
the scaling of a series of symmetric tensors, by different scale factors. This shows that the spin-s channels are never 
depolarizing channels. 

At the critical value p = A/3, the channel maps an arbitrary p v ^ w into a matrix with a Bloch representation: 

1 2 

£s(Pv,w)\ p=x/3 = ^l + -v-J. 

It follows that if p = A/3, then the channel maps an initial density matrix of the form (23) with v = into pure noise. 
We investigate critical values of p more systematically in Section V.F. 

For spin-1, an arbitrary density matrix may be represented as (23), and for higher spin, these are a proper sub- 
manifold of the convex cone of all density matrices. For spin 3/2, an arbitrary density matrix may be written in the 
form 

^ WabJ(aJb) + U abc J( a J b Jc) 
a, b a,b,c 

where w and u are completely symmetric tensors. The U-term is traceless, and so we require the W-term to have 
trace one. As discussed prior to Theorem 3, this means that tr(w) = 3/(dA) = 1/5. 



E. Finding v and w from p in Higher Spin and Pure States for Spin One 

In this brief subsection we show how to invert the relation (23) for the density operator, and find the coefficient 
vector v and symmetric tensor w. We do the analysis at arbitrary spin, although for spin higher than one, not all 
density matrices have the form (23). The methods will generalize assuming the relevant trace identities can be found. 

As in Theorem 5, let J\,J%, J3 be canonical generators for the spin-s representation of SU2. Note that 

Tr( J a J b ) = \d\S ab 

where A = s(s + 1) and d = 2s + 1. Also, 

Tr(J a J b J c ) = i ^ t abc 


where e a b c is the Levi-civita alternating symbol. It follows immediately from (23) that 

3 

v a = -TT-Tr(pJ Q ) . 
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To find w, note the trace identity 



Tr(J( a J b ) J(j Jfe)) = /l(s)-((5 a fe^ + 5 bk S a: j) + f2{s)S a bSjk ■ 



where /, are functions of s, given by 



/ 1 =tr(W= Md2 - 4) 



/ 2 = tr (Jfjl) 



30 

Ad(l + 2A) 
30 



By calculating tr(pj(j Jk)), we find 



30 / \ 2A+1 . xr 

For spin one, /i = 1/2, / 2 = 1, and rfA = 6 so 

1 1 

w a = j tr (P J a)' and = tr (p J U J k)) - 2 5 jk ■ 

This gives another way to find pure states: if p — \ip) (ip\ then 

I 1 

v a =^{Ja) i ,, and Wjk = {J(jJk))^-^jk- (30) 

The conclusion is that if the pure state p = \4>) has a representation of the form (23), then we can find its Bloch 
vector and w-matrix easily. 

Using the same conventions as in Section V.C for the spin one operators, equation (30) implies that for a pure state 
\ip) with three complex components ip a , we have 

Wab = \^ab - ^(tpalpb) ■ 

It also follows from (30) that 

v = ip R xip I , (31) 

where ipn denotes the real vector with components $l(ip a ), and ipi has components Im(^ ). The set of all v satisfying 
(31) with (ip | ip) = 1 is a ball of radius 1/2 in K 3 . 

It seems that there are no pure states for spin 3/2 in the space spanned by J a and J( a Jb)- To find the pure states, 
it is necessary to consider triples, i.e. u-states of the form 



abc 



F. Decomposition of the Space of Density Matrices into Convex Critical Regions 

In general, there is no value of p £ [0, 1] for which the spin-1 channel maps all initial density matrices into pure 
noise, so this channel is in some sense less decohering, and therefore more desirable, than its spin-half counterpart. 
Rather, there are two critical values, and two critical regions in the space T> of all density matrices. 

Note that we may re-write (23) in three dimensions as 

Pv,w = ^1 + V ■ J + ^ (w a b - ^5ob) J(aJb) ■ (32) 
a, b 
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The term containing the symmetric tensor is now traceless, and vanishes if and only if w = g 1. Define i?o to be the 
singleton set i?o = an( l 

i?i = {v ■ J : v £ M 3 } Pi = 2 

^2 = { Ea,h ^ ^ab) J (a^fc)} P2 = | 

In other words, elements of i?o> Ri, R2 are 3x3 matrices that take the respective forms of the three terms in (32). 
Note that R\ and R2 are convex sets containing only traceless matrices. 

Moreover, any 3x3 density matrix can be written uniquely as a sum of the form 

ro+n+ r 2 , where n £ R t . 

In compact notation, we have 

V(C 3 ) = R + R 1 +R 2l RiHRj = if i ^ 3, 

where pi is a critical value for R4, in the sense that any p £ Rq + Ri is mapped to pure noise at p = Pi- Density 
matrices not in Rq + Ri for some i will not have a critical value. 

This kind of decomposition holds for some (but not all) Lie algebras other than SU2. The spin-1 example (32) already 
well illustrates the fact that some of the critical values of p may lie outside the interval [0, 1] where the channel is 
trace-preserving. In that case, only p 2 = 2/3 is a true critical value, while p l — 2 does not describe a channel. In 
our general discussion of this phenomenon, which culminates in Theorem 6, we give an equivalent condition for the 
critical values to lie in the allowed interval. 

For any faithful, irreducible representation a of a semisimplc Lie algebra g on a Hilbert space H with d = dim(W), 
let £ denote the associated channel (5). Let X a denote generators of g, and also their images under a. We use the term 
rank to mean the degree of a polynomial in X a ; for example X a X b — ^\ ViX a has rank two. Rank is a well-defined 
function on the tensor algebra of g, but not on the universal enveloping algebra U(g), as Lie algebra identities relate 
polynomials of differing ranks. However, a given element of U(g) always has a (not necessarily unique) representative 
of minimal degree in the generators. 

Theorem 6. The space T>(H) of all density matrices on Ji admits a finite decomposition 

V{H) = R + Ri + --- + R N , 
K, lij if 1^ j, 

where Rq = {(l/d)l} ; each R r is a convex set consisting of traceless degree r combinations of the generators. Further, 
3 p r £ [0, 1] such that 

£(p) = ^1, at p = p r , for all p £ R n + R r 

if and only if the generators in this representation satisfy special r — > r — 2 identities with g r < 0. (It is most natural 
to define the term 'special identity' and to define g r below, following eqn. (38).) 

Proof. Let X a denote the representation matrices in the representation a. By surjectivity of a in the commutative 
diagram (18), we may write any density matrix p as 

p = v a X a + w ab X a X b + u abc X a X b X c + ... (33) 

a a.b a,b,c 

Let N be the smallest integer such that any p £ V(TL) can be written in the form (33) with at most terms. 
We may write any matrix T as the sum of its trace part and its trace-free part 

T = T tr + T 

where 



T te = - Tr(T) 1, and T = T - T tr . 
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By assumption, (33) has unit trace. It is then clear that the sum of the trace part of each term must equal (l/d)l. 
We may therefore rewrite (33) as 

a a.b 

Define i? to be the singleton set {^1}, and for n > 1, define R n to be the set of all matrices of the form 

W ai ,...,a n X ai . . .X an 



ai,...,a n 

,th 



for all w in the n symmetric power of R . It follows that the space of all density matrices is decomposed as 

V(H) = R q + Ri + --- + Rn, (34) 
Ri r\Rj = <D, i± j, 

Now let p r G Ro + Rr, so that 

Pr = -1 + ^2 S-fA-'-^-cl, (35) 
ai,...,a r 

for some coefficient tensor ui and constant c (equal to the trace of the rank r term). 

To simplify notation, we describe the relevant procedure for a rank 3 object p 3 G Ro + R3, with the understanding 
that the generalization to arbitrary rank is technically the same, but notationally worse. Writing 



P3 

(1 

abc 

where c is a constant chosen to make the terms in parentheses traceless, we then have 

k 



£(p 3 ) = (l-p)p 3 + ^Y.X*P3X*- (36) 



Recall that the representation-dependent constant Z is defined by the relation 



}Zxi = z-i 

i 

and is related to the quadratic Casimir and the Killing-norm of each of the generators. 
Suppose that the representation being studied has a 5 — > 3 identity, so that 

XiX a Xf,X c Xi = /o6cl + 9abcX a Xi,X c (37) 

for some tensors /, g. There is no implied sum on the rhs of (37). Consider using this to simplify (36), keeping only 
the degree 3 terms in the generators. The result is 



^ ^abc (l~P + ^gabc^j X a X b X c 



abc 

To make this expression vanish, we would like to solve the equation 

1 - p + | 5a6c = 0, (38) 

but that equation only yields a specific value for p when g abc = g 3 € C is a constant, i.e. takes the same numerical 
value for any selection of the indices a, b, c. 

This motivates the following definition of new terminology. In general, when an r — > r — 2 identity of the form (37) 
holds with g aoc equal to a scalar g r G C, let us call it a special r — > r — 2 identity. We have proven above that a 
special 4^2 identity (ie. g abc constant) exists for any irreducible representation of SU2, and found the form of that 
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identity. It is not hard to prove that for any irrep of SU2, special n — > n — 2 identities exist for all n. It is assumably 
an open question in representation theory whether they exist for other representations; we hope that the present work 
will motivate a further investigation of this important question. 

Assuming the special 5 — > 3 identity in the example of interest, we have 



+ 



^-i) P + i\Y,^ bc x a x b x c 



abc 

Note that u> and / are rank-3 tensors, and to ■ f denotes the full contraction Lo a bcf abc - The coefficient of the rank 3 
term vanishes at the value of p which sets the number in square brackets to zero. This critical value is 

Z 

P3 ■= 



Z-g 3 

This value of p is in the allowed range [0, 1] if and only if 173 < 0. The negativity of g sometimes holds and sometimes 
does not; for example, the 4— > 2 identity (29) would satisfy g < for dimension less than 3. In any case, this clarifies 
the point that a critical p r exists if and only if there is a special r — > r — 2 identity with g r < 0. 
Since at this value of p, the channel maps p 3 into d _1 l, we infer from (39) that 

u)- f = c(Z - g 3 ) . 

In the presence of a special r — > r — 2 identity, the channel maps Rq + R r to itself, i.e. 

£ (Rq + R r ) C Rq + Rr ■ 

This means that £ effectively looks like a depolarizing channel when restricted to Ro + R r - By choosing 

Z 



p = p r := 



Z - g r 



only the term proportional to the identity survives. Since the channel is trace-preserving, this term must be d 1 1, 
and we then have 

£(R + R r ) CR . 

□ 



G. Relation to the Werner- Holevo channel and a New Conjecture 

Datta (2004) has shown that the spin-1 channel at p = 1 is equivalent to the Werner-Holevo channel 

£{ P ) = ^- [ Hp)i- P t ) . (40) 

Recall that in our notation, Mq = \J\ — p 1, so taking p = 1 eliminates the identity from the set of Kraus operators. 
For p < 1 and for the spin s representation with s > 1, we may view the spin channel as a generalization of the WH 
channel. 

The Werner-Holevo channel became famous as a counterexample to the AHW conjecture (Amosov et al., 2000). 
We infer by Datta's equivalence that the spin-1 channel at p = 1 gives precisely the same counterexample to the 
AHW conjecture, stated below. Therefore, multiplicativity does not hold gcncrically in Lie algebra channels. Once 
it was established that the AHW conjecture docs not hold for all q > 1, it was natural to conjecture that it holds for 
1 < q < 2 (King and Ruskai, 2004), and this was recently proved for the WH channel by Alicki and Fannes (2004). 
If this is true, one would expect additional counterexamples with values of q approaching 2. However, none have yet 
been reported, except for the WH channel which gives a sequence of counterexample with q increasing from 4.79 as 
the dimension d increase. M. B. Ruskai, in a private communication to the author, suggested the possibility that Lie 
algebra channels might provide additional counterexamples with special properties: 
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Conjecture 1. Lie algebra channels generate counterexamples to the AHW conjecture (stated below) for a sequence 
of values of q approaching the boundary of the region in q-space where multiplicativity begins to hold for all channels, 
assuming there is such a region. 

For completeness, we now state the AHW conjecture, for which we need a definition. 

Definition 3. The maximal £ q -novm of a channel £ is defined as 

v q {£) = sup ||£( 7 )|| (g>l), 

where \\A\\ q = (Tr | A\i) 1/q , andV(H) denotes the space of density matrices on H. 

Amosov et al. (2000) conjectured that v q {£) is multiplicative for tensor product channels: 

v q {£® m )= sup ||£® m (r)|| = Vq {£) m . (41) 

Equation (41) is often called the t q multiplicativity relation or the AHW conjecture. Giovannetti et al. (2004) have 
conjectured that (41) holds for the Werner-Holevo channel when d > 2 q ~ 1 . 



VI. CHANNELS BASED ON EXCEPTIONAL LIE ALGEBRAS 

Let ej (j — . . . 7) denote the standard basis for the octonions O, where eo is the unit. Our notation is compatible 
with that of Baez (2002), and the proofs of our statements about the octonion algebra can be found there. The Lie 
group G 2 is the automorphism group of O, so the Lie algebra g 2 is the derivation algebra of the octonions: 

02 = 0et(Q). 

Derivations act trivially on the identity, and the imaginary octonions Im(O) form the fundamental 7-dimcnsional 
irreducible representation of q 2 . 

It is known that if A is an alternative, non-associative algebra (such as the octonions), any pair of elements x,y e A 
define a derivation D(x, y) : A — > A by 

D(x,y)a=[[x,y],a\-3[x,y,a\ (42) 

where [a, b, x] denotes the associator (ab)x — a(bx). When A is a normed division algebra, every derivation is a linear 
combination of derivations of this form. For the octonion algebra, the elements 

D(ei,e;), D(e 2 ,ej), and D(e 4 ,e k ) 

for all i > l,j > 2, and k > 4, are linearly independent and there are 14 such elements, so they are a basis for g 2 . 
Define the notation 

This is one possible basis for the Lie algebra g2, but we will use another more suited for our purposes. The fact (Mac- 
farlane, 2002) that the six-dimensional sphere S 6 may be viewed as a G2/SU(3) coset space, implies a corresponding 
decomposition of the algebra: 

£l2=m+J), f)=su 3 

where m is a 6-dimensional subspace. We find a basis adapted to this decomposition. The basis vectors for m are 
simply expressed as = di^+i, while 

hi = di2 + 2d 47 , h 2 = d 13 - 2d 46 , h 3 = d u - 2rf 2 7, 

hi = rfi5 + 2c?26, ^5 = ^16 — 2c?25, ^6 = ^17 + 2^24, 

h 7 = V3d 2 3, h s = d 2 3 + 2c?45 



are a basis for SU3. Let 



P = -j= ^{mi, ■ • • , me} U ^{^i- ■ ■ ■ > h s} 
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denote a corresponding basis for 02- Interestingly, [3 is an orthonormal basis of 02 with respect to the trace form on 
the 7-dimensional representation space, 

1 14 

Tri m(0 )(A/3j) = -Sij, therefore = h- 

i=l 

The 02 channel acts as 

14 
i=l 

Assume p has a Bloch representation with v € R 14 , 

p=i(/ 7 + «-/3), (43) 

then as an intermediate step, 

i=l 

The sum of /3? gives the identity, with a factor of p/7 to cancel the —p/7, and (miraculously) the term which is cubic 
in P's vanishes identically! This is due to the following remarkable 3^0 identity 

2&A.A = for all a, 

i 

as may be checked explicitly. Therefore, the 2 channel (restricted to its Bloch manifold) is the simplest of all. It is 
a true depolarizing channel, shrinking its Bloch vector by a factor of 1 — p, 



£{p) = \(l+{l-p)v-p) . 



We emphasize, however, that the 02 channel is almost certainly not a depolarizing channel outside the Bloch manifold, 
though we have not proven this. 

This does show that the critical value p\ = 1, in the notation of Theorem 6. 



VII. CHANNELS BASED ON THE CLIFFORD ALGEBRA 

Let ( , ) be a nondegenerate bilinear form on V, a <i-dimcnsional vector space. A representation of the Clifford 
algebra associated to (V, (, )) is a map 7 : V — > Ql(V) satisfying 

{l{x), 7(2/)} = (x,y)l. 

where the left side is an anticommutator. The representation is Hermitian if the image of 7 is contained in H(V), 
the (Hilbert) space of Hermitian operators on V. 

Theorem 7 (Clifford Algebra Channel). Given a Hermitian representation of the Clifford algebra, and a finite 
collection of nonzero vectors x\, X2, ■ ■ ■ , x n € R d , then 

n _ j n 

£ci{p) = (^(xuXi)^ ^2l{xi) pj{xi) (44) 

i=l i=l 

is a CPT map. 

Proof. The operator is completely positive because it is already in the form of an operator sum representation. Wc 
check that it is trace preserving. By cyclicity of the trace, 

n —\ n 

Tr(£ct(p))=(j2^ x j) £TKP7(z0 2 ) 

i=l i=l 

However, ~{{xi) 2 = ^{-f(xi),j(xi)} — (xi,Xi)l using the Clifford algebra. The sum of such terms decouples from 
Tr(p) and exactly cancels the prefactor. □ 



20 



We remark that, although the proof of Theorem 7 is trivial, the result may not be easily obtained by inspecting 
any of the standard matrix representations. Taking the Weyl representation of the 7 matrices in d = 4, one finds that 
writing out the CPT map 

- l (x)p 1 (x) + j(y)fry(y) 

for general x, y, p as an explicit matrix takes a full page. 

As we have seen in other examples, the computational methods used in this paper are most effective when an arbi- 
trary density matrix can be written in terms of the generators of the symmetry algebra. For the Weyl representation 
of the Clifford algebra, there is a convenient basis consisting of antisymmetric combinations of 7 matrices, which we 
summarize in the following table. 

1 1 of these 
7 M 4 of these 

7"" = ^[y,7"] = 7 [ 'V 1 6 of these 
y«v = 7 [a^ 7 p] = j e A^ 7ff7 5 4 of thcsc 
-r pa = ^yyy] = -ie^P ^ 1 of these 

These 16 matrices form a basis for the space g[(M 4 ). One can therefore write any 4x4 density matrix as a linear 
combination of these matrices with coefficients that are tensors of rank 4, and use 7 matrix identities to calculate the 
action of the CPT map (44). 



VIII. THE BLOCH MANIFOLD 

A. General Results 

The Bloch manifold is a geometrical space which is naturally associated to a certain representation of a semisimplc 
Lie algebra g, by asking the question: which linear combinations of the generators of g in that representation can be 
valid density matrices? For any preferred class of matrices (such as those with nonnegative eigenvalues) one can define 
a manifold from a representation in a similarly basis-dependent way, but for the application to quantum physics, we 
restrict interest to density matrices. 

Why is this an important question? We have already shown that the action of the su n channel is most simply 
expressed as a rescaling of the Bloch vector, and we know in that case that any n- level density matrix admits a Bloch 
representation. All that remains for a complete mathematical description of the su„ channel is to know the set of 
vectors v € R" _1 on which the transformation is being applied. 

What about other representations? The simplest example of why this is an important question for other represen- 
tations is the formula previously derived as (22), which gives the action of the spin-1 channel on density matrices 
admitting a Bloch representation (g = SU2) as: 

5( Pl) ) = I(l+(l-|) v . j) . 

Thus, the spin-1 channel also is a rescaling of the Bloch vector, and so it is not only for the n representation of su„ 
that characterization of the Bloch manifold (the space of admissible v) is important. 

In any Lie algebra representation which has a "triple product" or 3 — > 1 identity, i.e. an expression for 
J2i a(Xi)a(Xj)a(Xi) in terms of the generators a(Xk), any density matrix admitting a Bloch representation trans- 
forms very simply under the action of the Lie algebra channel. 

Definition 4. Choose a set of generators {Xi} of a semisimple Lie algebra g, and an irreducible representation 
a '■ ~* flK^O on a d a - dimensional Hilbert space H. The Bloch manifold V (in the Xi basis) is defined to be the set 
of vectors v G M fc (k — dimg) 7 such that 

p(v) = -^-(l + J2 v M X i))> ( 45 ) 

i 

is a valid density matrix. A density matrix which can be written in the form (45) is said to possess a Bloch represen- 
tation, and the corresponding vector v is said to be a valid Bloch vector. 
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In the notation of Theorem 6, the Bloch manifold is precisely the space Rq + R\ which appears in the natural 
decomposition of the space of all density matrices into critical regions. 

Theorem 8. The Bloch manifold is a closed set in R fc . 

Proof. The matrix p(v) is positive iff the lowest eigenvalue A min of 1 + J2 i Via(Xi) lies in the set [0, +oo). The lowest 
eigenvalue of a matrix is a continuous function of the matrix, so \ m in is a continuous function of v. The inverse image 
of the closed set [0, +oo) must be closed. □ 

Theorem 9. Let a be a d- dimensional representation of q, let k = dim(g), and let X a be an orthogonal basis of g 
with respect to the Killing form. Then 

Tr(a{X a )a(X b )) = Nd5 ab . 

Moreover, if v is a valid Bloch vector, then 

2 d - 1 

We remark that in the notation of Section III, 

N = Z/dim(g). 

Proof. The density matrix p = d _1 (l + v a a(X a j) must satisfy tr(p 2 ) < 1. But 

tr(p 2 ) = d-\l + Nv 2 ) (46) 
which implies the desired result. □ 

If {X'j} is a second basis of q, related to the original basis by a matrix A, then the Bloch manifold in the X 1 basis 
consists of A T applied to the Bloch manifold in the X basis. If det(j4) = 1, this yields an isometric copy of the 
original manifold, but otherwise the manifold has been stretched in some way. We will see examples of Lie algebra 
representations which are analogous to the qubit representation, in the sense that the Bloch manifold is a closed ball 
in some preferred basis. 



B. The Bloch Manifold for All SU2 Representations. 

As an example, we give the Bloch manifold relevant to the spin-j representation of SU2 • Let hj+i be the (2j + 1)- 
dimensional identity matrix, and the Ji are the standard generators in the spin j representation. The lowest eigenvalue 
of 

3 

^2j+l + Jj 
i=l 

is given by 1 — j \\v\\. We have proven: 

Theorem 10. The valid Bloch vectors for the spin-j representation 0/SU2 (with the standard basis) are elements of 
a closed ball in E 3 with radius 1 / j . 

Thus, the picture of the Bloch manifold as a closed ball is not necessarily particular to the qubit system, however, 
it is certainly not always a closed ball. As we shall see below, the Bloch manifold for the defining representation of 
su„ with n > 2 is a proper subset of the analogous closed ball. To complete the SU2 case, we remark that the radius 
receives a multiplicative constant if we rescale the generators; however, the radius always scales as one inverse power 
of the dimension of the representation. 
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C. A Bloch Submanifold from the Cartan Subalgebra 

The a(Xi) are Hermitian operators on H, which cannot in general be simultaneously diagonalized (if they can 
be, then cither g is abelian or the representation is trivial). Therefore, solving the positivity condition in more 
sophisticated examples is not straightforward. We discuss one method which works for any Lie algebra representation 
and which always gives a nonempty subset of the Bloch manifold. 

Let H\,...,H r denote a basis for the Cartan subalgebra of g, with r = rank(g). Simultaneously diagonalize all 
a(Hi), and let denote the j th diagonal element of a(Hi). The are, of course, weight vectors for the given 
representation. 

Assume that the basis {Xi} has the Cartan generators Hi, . . . , H r as its first r elements. We consider v £ R fe which 
are zero except for the first r components, which correspond to the Cartan generators, and ask when such a v gives 
rise to a density matrix. In this way we obtain a subset of the Bloch manifold. 

The positivity condition is 

r 

l + J2 V i h i j >° (Vj = l...da). 

i=l 

Each linear equation v ■ h 3 > — 1 defines a half-space H(j), and the restricted Bloch manifold 

V rcs = {v | v h? > -1 V j = 1, . . . ,d a } = f| M(j) 

is their intersection, clearly nonempty. For representations of nonabclian Lie algebras, 

V rcs C V. 

Depending on the rank of g, and on the spatial orientations of the weight vectors, the space V res is either a finite or 
a semi-infinite polyhedron. 



D. The Bloch Manifold for the Standard Rep of S\J n - 

We now discuss the structure of the Bloch manifold for the defining representation of su„. First, we derive a simple 
bound by applying Theorem 9 with N = 2/n, which yields 

9 n(n — 1) , , 

v 2 < -±-^-L . (47) 

By Descartes' rule of signs, an algebraic equation of degree n with real roots, 

n 

^{-iyajx n - j = \{{x - Xi ) = 0, x, g R 

j i=l 

has all roots nonnegative if and only if a* > for all i. It is then obvious that the Bloch manifold for the n-dimensional 
irrep of su„ is given by the set of v € K" _1 such that the characteristic polynomial ch p („\(a;) has only nonnegative 
coefficients. The coefficients a, can be calculated for any specific example using a = 1 and Newton's formula, 



a k = lj2(-l) q - 1 c q a k - g , (48) 



9=1 

where c q = J2 i x\ — tr (p(u) 9 ) . Naturally, calculating tr p q reduces to calculating the traces of products of at most q 
generators of su„. Since 

Oi = Ci = 1, and a 2 = ^(1 - c 2 ), 
the condition a 2 > is equivalent to trp 2 < 1, which leads to (47). 
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Using (48), the condition 03 > reduces to C3 > 5 (3c 2 — 1), but c 2 is given by (46), and a similar calculation shows 
that 

c 3 = trp 3 = (n + 6v 2 + 2v a v b v c d abc ) . 

The calculations up to this point have been valid for su„ for all n. However, to completely solve the problem for 
n > 3, we will need to know C4, C5, . . . Therefore, as a nontrivial example, we completely calculate the Bloch manifold 
for the SU3 channel in closed form. For n = 3, we note that 

2 

det(v ■ A) = -dijkVtVjVk 
so the condition c 3 > ^(3c 2 — 1) (for n — 3) can be expressed as 

det(uA) > -1 and v 2 < 1 + det(uA) 
Therefore, the Bloch manifold for the 3 of su 3 admits the following expression, beautiful in its simplicity: 

V SU3 = {v £ M 8 : v 2 < min (3, 1 + det(uA)) } . 

E. Bloch Manifold for the 7 of G2- 

In our calculation of the g 2 channel, we explicitly constructed a basis of g 2 using its definition as Det(O). This 
basis was normalized so that 

a 

Theorem 9 gives 

\v\ < 2V21. 

This proves that the g 2 Bloch manifold is contained in a closed ball of radius about 9.2. However, the true radius is 
much smaller, as we will now show. By g 2 symmetry, the 0's satisfy the identity 

tr (v0) q = 0, (Vd6 M 14 ), q odd 

where v0 = J2i v iPi- Further, for certain even values of q, tr (v0) q may have a simple expression. For example, 

tr (v0) 2 = V -, and (49) 
tr (v0f - (tr (v0) 2 ) 2 - ^ . 

These g 2 trace identity tr (w/3) 4 = (tr (v0) 2 ) 2 is not easy to prove. It is true because for g 2 and some other algebras, 
every fourth-order Casimir invariant is expressible in terms of the second-order invariant, as shown by Okubo (1979). 
Recently a simpler proof, together with other interesting trace identities, was given by Macfarlane and Pfeiffcr (2000, 
see (4.36)). 

Enforcing C3 > 7r(3c 2 — 1) gives a refinement, 

\v\ < 2V7 w 5.3. 

Requiring 04 > gives v 2 < 8(10 — \/65), so \v\ < 3.93. The coefficients are such that 05 > for all v, and tr (v0) n 
for n > 6 has no simple expression analogous to (49), so we have taken the simple analysis of the g 2 Bloch manifold 
as far as it will go. 
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F. Pure States in the Bloch Manifold 

Let the representation matrices be denoted by X a , a = 1, . . . , k. If v is in the Bloch manifold, so that 

Pv =d- 1 {l+J2 V aX a ) 

a 

is a density matrix, it is particularly easy to determine whether p is pure. If the products X a X b arc linearly 
independent from X a (i.e. there is no 2 — > 1 identity) then p 2 ^ p and the state is never pure. 

On the other hand, if the representation has a 2 — > 1 identity of the type satisfied by the fundamental representation 
of su„, 

X a X b = /3S ab l + ]T Q abc X c , (50) 

c 

then there can be pure states, and we have a complete characterization of them. 

Theorem 11 (Pure Bloch States). If the 2 — > 1 identity (50) holds, then a Bloch state p v is pure if and only if 

1 + flv 2 = d and ^ v a v b Q abc = (l - ^ v c , 

a. b 

for all c= 1, . . . , k. 
Proof. This follows from 

Pv 2 = ^2 I 1 + Z^ 2 ) 1 + ^2 J2( 2Vc + «a«6Qo5c)^c • □ 

abc 

It is interesting to see how Theorem 11 specializes to d = 2. In that case, Q a h c = and also 1 — 2/d = 0, so the 
second equation is always satisfied. The first equation amounts to v 2 — 1/(3, and f3 = 1, so this just says that w is on 
the boundary of the Bloch sphere, which is the well-known characterization of pure states. 

Unfortunately, 2 — > 1 identities almost never hold, excepting of course the fundamental representation of su„, 
because products X a X b tend to be linearly independent from the representation matrices X a if the dimension of the 
vector space is large enough to allow this. 



G. Summary of Bloch Manifold Technology 

The Bloch manifold in a certain basis is given by the solution of a system of polynomial inequalities in the com- 
ponents of the Bloch vector v. These inequalities come from enforcing positivity of the density matrix, p v > 0. In 
many cases, it is easy to see that the Bloch manifold is bounded within a ball, by enforcing the inequality tr(p 2 ) < 1. 
The Bloch manifold for the 3 of SU3 can be calculated exactly, and also in principle for $52- In the latter case, it is 
bounded within a ball of radius < 3.93. In any representation of any Lie algebra, if a 2 — > 1 identity (50) holds, then 
pure states lie on the surface of a sphere of squared radius (d — l)//3. 

What we have defined and studied here should rightly be called the linear Bloch manifold, because already for the 
spin-1 channel, one needs to represent the density matrix as v ■ J + J2 a b w abJ(a-h)- So the full geometry of the space 
of 3 x 3 density matrices is described by placing non-trivial conditions on both v and w, and similar remarks apply 
in higher dimensions. 

For the spin 3/2 channel, one describes the most general density matrix in terms of 

WabJ(aJb) + ^ U abc J^ a J b J c ) 
a,b a,b,c 

where w and u are completely symmetric tensors. Thus the full space of density matrices, in this representation, 
becomes a submanifold of the space of ordered pairs (w, u) G V® s2 ®V® s3 where V is the 4-dimensional representation 
space, satisfying some additional conditions. For a general representation V, to generate all density matrices one 
would need to consider a finite direct sum of symmetric tensor powers V®° n for various n. The answer becomes more 
complicated in higher dimensions (as does representation theory itself) but the method is completely general. 



25 



IX. CONCLUSIONS 

We hope that the reader will find it useful to have a compilation of results and formulae from the paper which lends 
itself to easy reference for future research. 

1. Definition of the generalized Lie algebra channel: 

k 

£{p) = (l-p)p + §£ a(Xi) pa{Xi) 

i=l 

where the basis is orthogonal, and Z is defined by X)i a (^) 2 — Z ■ 1. If Xi is Killing-orthonormal, then Z is 
the quadratic Casimir. 

2. Extension yields many other channels, including a "double g-channcl" with Kraus operators 



{^•V,.\\ : i,j = l...k}. 



These operators generate the image of g <g> fl under the universal homomorphism expressed in the commutative 
diagram (18). 

3. Action of the su„ channel on an arbitrary input density matrix: 



n \ n z — 1 



4. Minimal von Neumann output entropy S m i n of the su„ channel: 

np 



+ n 



l + n \n 2 -l) V 1 + nJ V 1 
with large n behavior: 

q . 

hm — — = p . 

n^oo m(nj 

This result was not discussed previously, but it is an easy calculation. 

5. An arbitrary 3x3 density matrix may be represented in the form 

Pv,w = v ■ J + S ^w ab J {a J b) (51) 

a.b 

for some vector v and symmetric tensor w, with tr(w) = 1/2, though not all objects of this form are density 
matrices. The action of the spin-1 channel on this density matrix is given by scaling the vector and tensor 
according to 

V a — > V a = (1 - -)V a , 
Wab -» W ab = (1 - —)W ab + -dab ■ 

6. Iteration of the spin-1 channel n times is the following transformation on w: 

Wab -p^ F (n) (p) (6 ab - 6w ab ) + w ab , 

where F^ n \p) is a degree n polynomial in p. F^(p) = 1 — 3p/2, and the for n > 1 are determined by the 
recursion relation 



^) W = (i-|)f(«)( p) + |. 
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7. Identities in the spin s representation of SU2: 

3 

J^JjJqJj = (A - 1) J a , where A = s(s + 1) 



i=i 

3 



^ JjJ(gJb)Ji — — 3) J (a Jb) + A5 a bl 



i=l 

Tr(J a J,,) = |dA5 o6 

Lr{J a J b J c ) = l — tabc 

^(J(aJb)J(]Jk)) = .fl(s)-^(5akSbj + S bk 8 a:j ) + f 2 (s)5ab8jk 



where fi are functions of s, given by 



. .2 Ad(d 2 -4) 
/1 = tr (J (1 J 2) ) = 
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8. Let Ji, J2, J3 be canonical generators for the spin-s representation of SU2 in dimension d — 2s + 1. For p„ )t(J = 
v ■ J + J2 a b w abJ( a Jb) to be a density matrix, we must have tr(w) = 3/ (cZA) where A = s(s + 1). If £ s represents 
the spin-s channel, then 



£ s (Pv,w) = (l - j)v ■ J+ (l - y)w ao J (a J fc) + ptr(w)l . 



9. With d, A, and p VjW as above, we have 

3 

v a = -j^tv(p v , w J a ) , and 

30 i t j \ 2A+1 

10. If a density matrix p w is pure and takes the form 

Pw — ^ w a bS( a S b ), 

a,b 

(i.e. v = 0) then the matrix elements are (p w )ij = icacij, where 

ai = V^22 + w 33 , a 2 = \Zj - w 22, «3 = \J\-wzd, ■ 

11. For the pure state V with components Va, we have 

w a6 = i<5 ah - $t(il>ai>b), and = 

where Vfl denotes the real vector with components ^t(ip a ), and V/ f° r the imaginary part. The set of all v 
satisfying (31) with (tp \ ip) — 1 is a ball of radius 1/2 in R 3 . 

12. The space V(7i) of all density matrices on H admits a finite decomposition 

V(H) = R + Ri + --- + Rn, 

Ri n R 3 ■. = if j, 

where Ro = {{l/d)l}, each R r is a convex set consisting of traceless degree r combinations of the generators. 
Further, 3p r € [0, 1] such that 

£{p) = ^1, at p = p r , for all p e i? + -Rr 

if and only if the generators in this representation satisfy special r — > r — 2 identities with g r < 0. Special 
identities were defined following eqn. (38). 
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13. Action of the g 2 channel on a Bloch- vector input: 

14. When <f>, defined by the commutative diagram 

^W(fl) 




is surjective, then the calculational technique used in this paper will always work. This surjectivity holds under 
a very general set of assumptions. 

15. If 7 denotes a representation of the Clifford algebra associated to the bilinear form ( , ) then the following 
expression defines a channel: 

i=l i=l 

16. The Bloch manifold is defined to be the set of vectors v e K fe such that 

P (v) = ^-(l + ^vMXi)) ev, 

i 

It is always a nonempty closed set in R fc , k = dimg. Finding this manifold is important because it parameterizes 
the space of density matrices for which the Lie algebra channel has a simple formula. 

17. In a d-dimensional representation of 0, normalized so that Tr (a(X a )a(Xb)) = dNS a b, for v in the Bloch manifold 
we have 

2^-1 

V < . 

- N 

18. In a certain basis, the Bloch manifold always contains the intersection of the half-spaces v ■ h? > — 1, where h? 
are the weight vectors of the representation. 

19. The valid Bloch vectors for the spin-j representation of SU2 (with the standard basis) are elements of a closed 
ball in R 3 with radius 

20. The Bloch manifold for the 3 of SU3 admits the following beautiful expression: 

V SU3 = {v € K 8 : v 2 < min (3, 1 + det(uA)) } . 

21. Vectors v in the Bloch manifold for the Q2 channel, with the chosen normalizations, satisfy 

v 2 < 8(10- V65), or \v\ < 3.93. 

22. If the 2 — > 1 identity (50) holds, then p is pure if and only if 

l + [3v 2 =d and ^v a v b Q abc = {l -^jv c , 

a, b 

for all c = 1, . . . , k. 
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